Interrelation between the distributions of kinetic energy release and emitted electron 

energy following a decay of electronic states 
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In an electronic decay process followed by fragmentation the kinetic energy release and electron 
spectra can be measured. Classically they are the mirror image of each other, a fact which is 
often used in practice. Quantum expressions are derived for both spectra and analyzed. It is 
demonstrated that these spectra carry complementary quantum information and are related to the 
nuclear dynamics in different participating electronic states. Illustrative examples show that the 
classical picture of mirror image can break down and shed light on the underlying physics. 
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The fragmentation of molecules and clusters following 
electronic excitation is a broad and widely studied pro- 
cess [ij-Q . The detailed knowledge of the process includ- 
ing the distribution of the produced fragments is relevant 
for atmospheric chemistry, for astrophysics and for basic 
research. The electronic excitation by photon or electron 
impact leads to excited electronic states which can disso- 
ciate by themselves or decay by emitting a photon or an 
electron, and the final states of this decay may then frag- 
ment. Photodissociation Q provides a prominent class 
of examples of the former case, and the fragmentation of 
molecules which have undergone Auger decay [S] forms 
an important class of examples of the latter. 
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We shall concentrate here on the fragmentation fol- 
lowing the decay of the excited electronic state and, to 
be specific, consider decays by emission of an electron, 
although, the situation is similar when a photon is emit- 
ted. Assume for instance a diatomic molecule AB where 
a core electron of A is removed by a high energy photon. 
The resulting A^B will undergo Auger decay and emit 
an electron, the Auger electron. If dications A^B^ are 
produced, they will, of course, be subject to a fi'agmen- 
tation into A^ and B^ via Coulomb explosion [gl, [J. An- 
other interesting whole class of processes is interatomic 
Coulombic decay (ICD) 0415 1. Assume an atom M and 
a rather distant neighboring atom N which is not chem- 
ically bound to M, an example would be a noble gas 
dimer. After removing a core electron from M, the Auger 
decay is essentially atomic and M^^ is produced. ICD 
then takes place as a follow up of Auger resulting in 
the - N+ trmle ion and an emitted electron, the 
ICD electron [11, [H-il | . The ion, of course, undergoes 
Coulomb explosion and fragments. 

In such experiments one measures the distributions of 
the total kinetic energy released by the fragments (KER 
spectrum) and the kinetic energy of the electron emitted 
in the decay (electron spectrum) which can be obtained 
separately or from the coincidence measurements of all 
charged particles [III, 113, HI il7l - Ell |. The question im- 
mediately arises whether at all the KER and electron 
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FIG. 1: (Color online) Potential energy curves of a system 
which undergoes an electronic decay process with subsequent 
fragmentation. In a classical picture the system decaying at 
time t at an internuclear distance i?o will produce emitted 
electrons with kinetic energy Ee{Ro) and the fragments with 
KER energy EKEn{Ro) = Et — Ee{Ro). Two sources con- 
tribute to the KER: the energy released from the potential dif- 
ference AVf — Vf{Ro) — Vf{oo), and the kinetic energy i?kin. 
In the quantum picture the wavepacket tpd propagates and 
continuously decays onto the final potential curve V/ where 
the resulting wavepacket tpf propagates. The KER and elec- 
tron spectra are then given by Eqs. (|4l5p . 



spectra contain complementary information on the sys- 
tem and the underlying process. Let us start discussing 
the problem in the framework of classical mechanics. Let 
the excited molecule possess the total energy Et' which 
consists of the kinetic energy i?kin of the relative nuclear 
motion and the potential energy Vd- At a specific time 
t the excited molecule decays at an internuclear distance 
Rq and emits an electron of energy Ee{Ro) which is given 
by the energy difference Vd{Ro) — Vf{RQ) between the po- 
tential energy curves Vd of the decaying and Vf of the fi- 
nal electronic states. The situation is depicted in Fig. [T] 
The produced ion now moves on the dissociative curve 
Vf and fragments. The KER energy Eker is obviously 
the difference between the energies of the ion produced 
at Rq and that of the fragments after they have com- 
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pletely separated: -Bker = [-E^kin + Vf{Ro)] - Vf{oo). 
The consequences are clear. If we know the potential en- 
ergy curves Vd{R) and Vf{R) of the decaying and final 
electronic states, there is nothing new in the KER spec- 
trum crKER(-E'KER) Compared to the electron spectrum 
ae{Ee) and vice versa. Because of energy conservation 
at every point R, one spectrum can be expressed as the 
mirror image of the other: ctkerC-E-ker) = cr^iET ~ E,,) , 
where Et is the total energy available relative to the to- 
tal electronic energy of the separated atomic fragments, 
i.e., Et = -Eker + Ee = Et' - Vf{oo). 

Indeed, the mirror image procedure has been success- 
fully verified in experiments and also computationally, 
11 , 23] • Are there no fingerprints of quantum 
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dynamics which make the spectra inequivalent? We shall 
show below that quantum dynamics can easily lead to a 
different physical content of the two spectra and then the 
mirror image procedure breaks down. 

We proceed by deriving quantum expressions for the 
KER and electron spectra starting from the coincidence 
spectrum cr(£'KER7 E^, t) where all particles are measured 
in coincidence as a function of time. The KER spectrum 
is determined from the coincidence spectrum by integrat- 
ing over all electrons 

'T'KErI-E'kER, i) = y d£^e Cr(£^KER, -Ee, 0- (1) 

and the electron spectrum by integrating over all -Eker 



(^e{Ee,t) ~ j dEKER Cr(EKER , Eg , . 



(2) 



In current experiments the spectra are measured after a 
long time, i.e., t — )■ oo, and we will in our explicit ex- 
amples concentrate on this standard case. We would like 
to point out, however, that time-dependent spectra carry 
much more information on the physics of the underlying 
processes, see e.g. {24l |. 

One can further evaluate Eqs. (|ll2p by noting that 
the coincidence spectrum is determined by the nuclear 
wavepacket ipf{Eejt) propagating on the potential en- 
ergy curve Vf{R) of the final electronic state to which 
the system has decayed [ll] . Before the decay has started 
this wavepacket is, of course, unpopulated, i.e., its value 
is zero. As time proceeds, tpf{Ee,t) becomes populated 
and carries all the information of the decayed system. Its 
equation of motion follows from the Schrodinger equation 
and is well known :25] : 

i\i;fiE,,t))^Wd^f\Mt)) + iHf+E,)\^bfiE,,t)). (3) 

Hf is the common nuclear Hamiltonian governing the 
nuclear motion on the potential curve V/(i?), and the 
transition matrix element from the decaying state to the 
final state is denoted by Wd^f ■ The initial excitation of 
the decaying state is, of course, provided by experiment. 



Once excited, the nuclear wavepacket ipd{t) of this elec- 
tronic state propagates on the complex potential curve 
Vd — iT/2, where F is the total decay rate of the state 
and may depend on R. With time electrons are emit- 
ted and ipdit) decays and loses its norm thus populating 
^pf{Ee,t) as can be seen in Eq. ([3]), in which Wd^fipd(t) 
is the source term for populating the final state. There 
is an intimate relation between F and Wd~^f- If there 
are several final electronic states / into which the system 
decays, then F/ = 2n\Wd^f\'^ is the partial decay rate 
and the total rate is simply F = ^^F/. In the follow- 
ing we assume for simplicity of presentation a single final 
state. We note, however, that the extension to several 
final states is trivial. 

The coefficients in the expansion of ipf{Ee,t) in the 
complete set of dissociating nuclear wavefunctions \Ef) 
of the potential Vj (R) determine as usual the probability 
of finding the fragments with energy Ef and the electron 
with energy Eg and thus determine the coincidence spec- 
Explicitly, \^pfiE„t)) = j:cEf{Ee,t)\Ef), 
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trum 

where \Ef) is an eigenfunction oi Hf with energy Ef — 
Eker + immediately leads to cr(£'KER, -Ee, — 

\cEj{Ee^t)\^ . Inserting this expansion into Eq. ^ and 
projecting on a single energy allows one to find the so- 
lution of the expansion coefficient for that energy. Inte- 
grating over the electron energy E^ gives via Eq. ((U an 
interesting expression for the KER spectrum: 



CTKER(-BKER,i) = 2^ / dt'\{Ef\Wd^f\i^d{t'))\ 



(4) 



This general expression is easily interpreted. First, the 
KER spectrum is solely determined by the dynamics in 
the decaying state and not at all by the dynamics in the 
final state. Second, as in any transition, the integrand 
can be viewed as a generalized Franck-Condon factor con- 
necting '0d and the dissociative eigenfunction of the final 
state. Third, and importantly, since V'd is time depen- 
dent, the KER spectrum at time t is given by the Franck- 
Condon factor which has accumulated up to this time, or 
briefly, by the accumulated Franck-Condon factor. 

Similarly, we obtain for the electron spectrum 



ae{E,,t) = 01^-/(^^6,0), 



(5) 



which is, however, a well-known result |25l. l27l|. In sharp 
contrast to the KER spectrum, the electron spectrum is 
determined by the norm of the final wavepacket which 
depends on the energy of the emitted electron, and thus 
solely refiects the dynamics in the final state. Of course, 
this dynamics itself is not independent of that in the de- 
caying state as can be seen in Eq. ([3]). Accordingly, ij^dit) 
evolves and decays and its losses Wd^f'pdit) are contin- 
uously transferred to the final state where they continue 
to propagate but now with the final state Hamiltonian 
Hf . These transferred losses can, in principle, collide 
with the already available parts of ipf{t) and interfere 
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with them. All of these rather complex dynamical phe- 
nomena are missing in the accumulated Franck-Condon 
factor which determines the KER spectrum. In this re- 
spect the electron spectrum carries much more intricate 
information. We would like to remind, however, that the 
differences between the two kinds of spectra are due to 
quantum effects. 

After having derived and interpreted the full quantum 
expressions for the KER and electron spectra we present 
illustrative examples. A transparent model example and 
two realistic examples which have been recently mea- 
sured. In the model example depicted in Fig. [U the de- 
caying state is described by a harmonic potential curve 
Vd (vibrational frequency 200 meV) and the final state of 
the decay by the dissociative curve 2/R + Vf (oo) which 
is chosen to be the same as in one of the realistic ex- 
amples discussed below. The decaying electronic state 
decays by a constant rate F = 200 meV which is within 
the typical range of Auger decay rates [H, [5^ . The sys- 
tem is initially in its ground electronic (and vibrational) 
state which is also chosen to have a harmonic potential 
curve but which can be shifted with respect to Vd- As 
often done, the system is excited by a broad band light 
pulse transferring its initial vibrational Gaussian wave- 
function V'd(O) vertically to Vd- This wavepacket ipd{t) 
can swing to and fro on Vd and thereby decay as elec- 
trons are continuously emitted. The parts of ipd{t) lost by 
the decay {d{'ipdit)\ipd{t)) /dt = -F(V'd(t)|'0d(i))) appear 
as ipf{Ee^t) on the final dissociative curve, see Eq. ^ 
describing the dissociation dynamics of the system. 

The results of the numerical calculations on the model 
example are collected in Fig. [2] Let us begin by choos- 
ing equilibrium distances of the ground and decaying 
states to be the same. Then, only the lowest vibra- 
tional level of Vd is populated by the optical transition. 
Correspondingly, no nuclear dynamics takes place on Vd, 
and one has V'd(0 = '0d(O)e~^*/^. What to expect in 
such a simple situation? Eq. (jlj can be easily solved 
giving ctkerC-EkeRjOo) = |(-Eker + V/(oo)|?/'d(0))|^, i.e., 
the KER spectrum is determined by the usual Franck- 
Condon factors. The decay rate F does not influence 
at all the KER spectrum! On the other hand, even in 
this simple situation, this is not the case for the electron 
spectrum, which turns out to be given by convoluting 
the mirror image of (Tker(-E'ker, oo) with a Lorentzian 
of full width at half maximum (FWHM) F. The nu- 
merical results are shown in Fig. [2^ together with the 
KER spectrum obtained as the mirror image of the elec- 
tron spectrum as expected classically (denoted 'predicted 
KER'). Clearly, the exact KER spectrum is narrower and 
exhibits a nodal structure at its high-energy wing being 
a fingerprint of the dissociative eigenfunctions of state /. 
According to the above, convoluting this KER spectrum 
with a Lorentzian of FWHM F reproduces exactly the 
classically predicted KER. 

Shifting the equilibrium distance of the system to 
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FIG. 2: (Color online) The electron and KER spectra of a 
model study (see text and Fig. 1). The initial wavepacket is 
chosen to be a Gaussian (lowest vibrational level of the ground 
electronic state of the system) in both panels: in panel (a) it 
centers at 3.5 A, in panel (b) at 3.3 A. F and the potential 
curves are given in the text. These spectra are computed via 
Eqs. (|4l5p in the limit f — >■ oo. The results are compared 
with the (classically) predicted KER spectrum obtained by 
mirroring the electron spectrum. The differences between the 
predicted and exact KER spectra (shown at the left-hand side 
in both panels by the black dotted and red lines, respectively) 
demonstrate that quantum effects can be substantial. 



smaller values ( 3.3 A in Fig. [TJ, introduces nuclear dy- 
namics of ifid on the decaying curve Vd- The decay pro- 
cess becomes more intricate even for a constant F and 
harmonic curves, and the electron spectrum exhibits vi- 
brational interference effects which have been discussed 
theoretically and measured [soj. The computed electron 
spectrum in Fig. O) indeed shows a typical vibrational 
progression including the impact of the interference ef- 
fects. Interestingly, the exact KER spectrum possesses 
a very different structure and does not at all resemble 
the classically predicted picture of mirror imaging the 
electron spectrum. Instead of a progression of declining 
peaks ending with a long low-energy tail, the exact KER 
spectrum essentially consists of a small and a pronounced 
sharp peak at low energy and a long high-energy tail. 
Different from the electron spectrum, where the struc- 
ture shows the fingerprint of different vibrational levels 
of Vd, the sharp peak together with the small shoulder of 
the KER spectrum are contributed by ^d at the classical 
turning points. This is verified by evaluating the time 
evolution of Eq. Q with a semiclassical approximation 
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FIG. 3: (Color online) Computed and experimental KER 
and electron spectra for two different processes. Upper pan- 
els: Auger process in CO (C0+ C0'^+{2^j:+) + e ). Lower 
panels: ICD in NeAr following the Auger decay of Ne^ Is"'^ 
(Ne2+(2s-i2p-i ^P)Ar ^ Ne'+{2p-' ^D) - Ar+(3p-i) + e"). 
(a) Auger KER spectra. The experimental spectrum is from 
Ref. [3j|. To compare with experiment, the computed KER 
spectrum (red curve) is convoluted with a Gaussian of 0.15 
eV FWHM. (b) Auger electron spectra. The experimental 
spectrum is from Ref. [s^]- Shown are also two computed 
spectra. The spectrum as calculated using Eq. [S] (blue curve) 
and that obtained by convoluting it with a Gaussian of 0.68 
eV FWHM to account for the experimental resolution (or- 
ange curve). Clearly, the electron and KER spectra are far 
from being the mirror image of each other, (c) Computed 
ICD spectra. Shown are the KER spectrum calculated via 
Eq. |4] (red curve) and the ICD electron spectrum computed 
using Eq. [S] and mirror imaged to provide a KER spectrum 
as predicted by the classical theory (black curve), (d) The 
experimental KER spectrum [35l ] is compared with the com- 
puted KER spectrum obtained by convoluting the black curve 
in panel (c) with a Gaussian of 0.7 eV FWHM [ssl] to account 
for the experimental resolution (black curve). Theory and 
experiment compare well. In this example the prediction of 
classical theory holds that the electron and KER spectra are 
the mirror image of each other (see also text). 



for \Ef). 

Let us now examine realistic examples, namely the 
Auger decay of CO and ICD of NeAr. Consider the Auger 
channel CO+(C Is'^) ^ C0^+(2 + e", whose elec- 
tron and KER spectra have been measured SlJ, 13211 . Po- 
tential curves and F are taken from literature 1331 13411. 



The experimental and computed KER spectra depicted 
in Fig. [Sja) coincide well. The computed spectrum is 
computed via Eq. ^ and broadened to account for the 
experimental resolution. This procedure is in accord with 
the discussion of the fragmentation process in 
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sure this spectrum, the resolution is much lower than for 
the KER spectrum [3l|, Isl] and we show both the spec- 
tra computed via Eq. ^ (blue curve) and the one convo- 
luted with a Gaussian resolution function (orange curve). 
The latter reproduces the experiment. Regardless of the 
experimental resolution, it is obvious from these calcula- 
tions that the KER and Auger electron spectra are far 
from being the mirror image of each other as expected 
classically. 

Contrary to the molecular Auger process, the electron 
and KER spectra obtained for the ICD processes in noble 
gas dimers are usually considered to be the mirror im- 
ages of each other. Here we study as an example of cur- 
rent interest the ICD process Ne^"'"(2s~^2p~^ ip)Ar 
Ne^+(2p-2iD) - Ar+(3p"^) + e~ following the Auger 
decay of Ne'''(ls~^). This very fast Auger decay pro- 
duces Ne^'*' atomic dicationic states and some of them 
can further decay by ICD producing Ne^'^'Ar"'' triply ion- 
ized states which undergo Coulomb explosion. The pro- 
cess has recently been measured 1351. In our calculations 
we utilize the potential curves of [22| and the value of F 
suggested in [35| . The mirror image of the computed ICD 
electron spectrum and the computed KER spectrum are 
shown in Fig. ^c) and agree well with each other. This 
agreement is due to the facts that the potential Vd is very 
shallow and the rate F is rather small. Consequently, sev- 
eral quasi-degenerate vibrational levels of Vd are initially 
populated, and the system behaves as if a single effec- 
tive level is populated. A comparison with experiment is 
provided in Fig. ^d). 

Although being the mirror image of each other within 
the classical picture, the KER and electron spectra of 
a decaying state carry complementary quantum infor- 
mation on the decay process. While the latter reflects 
the nuclear dynamics in the final electronic state and is 
sensitive to interference effects, the former measures the 
accumulated Franck- Condon factor of the decay and is 
the projection of the dynamics in the decaying electronic 
state on the potential curve of the final state. The ex- 
plicit general expressions derived allow to compute and 
analyze these spectra, especially if a pulse is involved. 
Illustrative examples show that the classical picture of 
mirror image can break down due to quantum effects. 
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The experimental and computed Auger electron spectra 
are shown in Fig.[3Ub). Because of the difficulty to mea- 
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